%有种群竞争 60分钟后湿度和温度会动
function dx = func32(t, x, wa, T) 
    global w v 
    dx = zeros(3, 1); 
    H1 = 0.8; H2 = 0.4; H3 = 0.2; 
    N1 = 2000; N2 = 2000; N3 = 2000; 
    sigma21 = 0.7; sigma32 = 0.7; sigma31 = 0.56; 
    sigma12 = 0.8; sigma23 = 0.8; sigma13 = 0.64;    
    k1 = 2; k2 =2.2; 
    
    %if t>60
    %    T=32+5*sin(0.1*t);
    %    wa=0-2.5*sin(0.2*t);
   % end
%T(ceil(t))
    u = water(wa); y = temp(T(ceil(t))); 
    %y=temp(T);
    sigma = k1*u*(1-w)+k2*y*(1-v); 
    dx(1) = H1*x(1)*(1 - x(1)/N1 - sigma21*x(2)/N2 - sigma31*x(3)/N3 - sigma(1)); 
    dx(2) = H2*x(2)*(1 - x(2)/N2 - sigma12*x(1)/N1 - sigma32*x(3)/N3 - sigma(2)); 
    dx(3) = H3*x(3)*(1 - x(3)/N3 - sigma13*x(1)/N1 - sigma23*x(2)/N2 - sigma(3));  
end